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Abstract 



We present a theoretical study for the intermediate stages of the growth of 
' membranes and vesicles in supersaturated solutions of amphiphilic molecules. 

The problem presents important differences with the growth of droplets in 



the classical theory of Lifshitz-Slyozov- Wagner, because the aggregates are 
extensive only in two dimensions, but still grow in a three dimensional bath. 
\^ I The balance between curvature and edge energy favours the nucleation of 

' small planar membranes, but as they grow beyond a critical size they close 

' themselves to form vesicles. We obtain a system of coupled equations describ- 

^ . ing the growth of planar membranes and vesicles, which is solved numerically 

. for different initial conditions. Finally, the range of parameters relevant in 

a 
o 
o 



' experimental situations is discussed. 

o ■ 



I. INTRODUCTION 

Solutions of amphiphilic molecules in water may form a large variety of molecular aggre- 
gates, as a result of the asymmetric interaction with water of the hydrophilic heads and the 
hydrophobic tails of these molecules . These aggregates range from micelles, which may 
be considered as large molecular clusters, to continuous structures of the amphiphile with 
regular or irregular structures such as sponge phases. In between, one finds structures like 
membranes and vesicles which are macroscopic in two dimensions, whereas they are only a 
few molecular sizes thick. From a thermodynamic point of view one can regard membranes as 
two dimensional phases in coexistence with a diluted bulk solution of water and amphiphile 
molecules, but contrary to the case of adsorbed layers on solid substrate, the membrane is 
not restricted to lie on fixed positions and orientations. The membrane is free to explore 
the full three dimensional space, and in the case of "fluid membranes", the lack of rigid 
molecular order within the membrane surface allows the bending of the membranes, with 
low energetic cost as "curvature energy" . The extraordinary properties of membranes as 
self-assembling, self-sealing, insulating and flexible structures are used in biological systems 
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to form the basis of cellular membranes, plasts and mitocondrial structures, the insulating 
material for neurons, etc. In the last decade, the study of membranes and other amphiphilic 
aggregates has attracted the attention of experimental and theoretical physicists, chemists, 
physical-chemists, biologists |7HT0|. 



In the present paper we discuss the dynamics of formation of these aggregates from a 
supersaturated bulk solution, by extending to the peculiarities of these systems the methods 
developed for "regular" three dimensional phases. We focus our interest on amphiphilic 
systems forming two types of isolated bilayer aggregates: membranes and vesicles. During 
the formation of these macroscopic aggregates we may distinguish two different stages: 
(i) nucleation process, during which the microscopic aggregates form. 
(a) A growth process, when these seeds reach mesoscopic sizes or dissolve into the bulk so- 
lution, by following a nearly deterministic dynamics which transforms the initial population 
of the microscopic clusters into a distribution of sizes for the aggregates and depletes the 
amphiphile concentration in bulk solution to the equilibrium value. 

In an infinite system at zero temperature, the growth process would never end. Alterna- 
tively, we should consider a third regime when the super-saturation is so small that the 
gradients in the chemical potential are affected by fluctuations. 

Experimentally, these aggregates can be obtained with the help of techniques designed 
to accelerate the nucleation process. Among these are the ultrasound technique and the 
use of solid substrates (heterogeneous nucleation), which preferentially adsorb the am- 
phiphilic molecules. On the other hand, the spontaneous formation from an equilibrium 
super-saturated solution is of interest. In fact, it has been suggested [|12| that, due to the 



possibility of spontaneous formation, amphiphilic systems forming vesicles played a crucial 
role in the beginning of life, providing for the compartmentation of the early metabolic and 
self-reproducing molecular machinery. 

In this work we address the growth process, in which the aggregates formed in the early 
stage compete with each other to grow, incorporating the amphiphilic molecules from the 
super-saturated solution. We try to answer the questions about how the typical distribution 
of sizes evolves and how it depends on the distribution in the early stage. Our model 
needs as input the initial distribution of sizes as well as the initial super-saturation of the 
system. This data could be obtained from nucleation theory, but the study of the different 
mechanisms of nucleation requires molecular models of these complex systems, which are 
still in an early stage of development. 

II. THE MODEL 

Let us consider the possible equilibrium shapes for bilayer aggregates. For a fixed area, 
A, a symmetric bilayer adopts the shape that minimizes the effective Helfrich Hamiltonian 

n = JdA [l^^{Cl + + KGCiC2j + Xfdl, (1) 

where the first integral extends over the whole surface area and the other one extends 
along the boundary, where n and are respectively the bending rigidity and the Gaussian 
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curvature moduli, A is the line tension and Ci, C2 are the two principal local curvatures of the 
surface. In the case of a flat circular membrane of radius R, Ti. oc R, whereas for a spherical 
vesicle of radius S, H. = 47r(2K + Kg). Thus, small aggregates tend to assume a planar shape 
but larger aggregates will always prefer a closed shape. To locate the spontaneous shape 
transformation we employ a simple spherical cap model where the bilayer is restricted to 
adopt the shape of the curved surface of a sphere cut by a plane. For a fixed area the shape 
can be fully described by the radius, R, of the circle of intersection between the sphere and 
the plane. 



where = 2k + hq- A spherical vesicle corresponds to i? = while a planar membrane is 
described hjR = (A/vr)^/^. In Figure 1 we display the energy H versus R for different values 
of A. The absolute minimum is always at an extremum: for A < AirlKs/Xy it corresponds 
to planar membranes while for A > 47r(/ts/A)^ it corresponds to spherical vesicles. When 
the two configurations have the same energy {A = 47r(Ks/A)^), they are separated by an 
energy barrier ATi = Sttk^, which decreases with increasing A and finally disappears at 
Ac = 7r(4fi;s/A)^, where all planar membranes become unstable and close into vesicles. 

Such a simple model gives a first glimpse of the growth process: unless the nucleation 
mechanism is strongly biased to produce vesicles, the early population of aggregates consist 
of planar membranes; as some of these grow they become metastable with respect to the 
closed vesicles, but still separated from them by a barrier, until their area reaches the critical 
value Ac- At this stage the membranes become unstable and turn into closed vesicles of the 
same area, which can grow until the solution has been depleted down to the equilibrium 
concentration of amphiphiles. 

We start constructing our model of growth for a system of both vesicles and membranes 
by making the following assumptions: 

(i) The excess concentration (with respect to a dilute solution without aggregates) is so 
small that the interaction between aggregates can be neglected. 

(a) All aggregates are much larger than their typical thickness ~ 40A). This assumption 
is expected to hold at least in the long time regime. 

(in) The main mechanism for transport of amphiphiles is diffusion as expected when the 
water solution is at rest, so that all hydrodynamic effects can be neglected. In this sense it 
is implied that the growth process does not agitate the water solution. 

(iv) We consider an infinite volume. 

(v) Finally we assume that fluctuations can be neglected. 

We shall later comment on some of these assumptions. Within this approximations the 
problem clearly resembles the classical Lifshitz-Slyozov- Wagner (LSW) theory for the growth 
of spherical domains 

Within a quasi-equilibrium thermodynamic description the chemical potential, /i(r), 
changes smoothly with position. We consider an isolated growing aggregate of area A and 
N{A) molecules; its rate of growth will depend on the number of amphiphilic molecules that 
approach its surface: 




(2) 




(3) 
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where the integral extends over a surface enclosing the aggregate and j(r) is the current of 
amphiphilic molecules induced by the gradients of chemical potential which is different in 
the bulk and at the surface of the aggregate: 

j = -aV/i(r), (4) 

a being a kinetic coefficient. Assuming that the transport of amphiphiles does not change 
its local concentration, 

V-j = 0, (5) 

everywhere outside the border of the aggregate. This leads to a Poisson equation for /i(r), 
equivalent to an electrostatic potential with boundary conditions at infinity, /^oo, and at the 
surface of the aggregate which mimics an equipotential metallic boundary. Such a 

value, depends on the particular properties of each aggregate. 



A. Membranes. 

Let us consider a system composed only of membranes and assume that none of them 
transforms into vesicles. The excess of grand potential energy of a planar circular membrane 
of radius R is 

^£lm[R) = TiaR^ + 27r\R, (6) 

where a is the surface tension of the membrane and A is the line tension associated with the 
boundaries. The equilibrium of a large membrane of arbitrary shape is controlled by the 
requirement a = 0, which at a given temperature is satisfied when the chemical potential 
/i assumes its equilibrium value /ig. In general, for a fixed chemical potential /x, the surface 
tension does not vanish but is given by: 

^(/^) = -r(/^ -/^o) (7) 

where F is the adsorption per unit area in the membrane. Since we are considering systems 
with small supersaturations, within a linear approximation we take the values of A and F 
evaluated at /i = /Xo. For /i > /io eq. @ has a maximum at 

In other words, membranes with R > RdfJ^) tend to adsorb particles and grow. This growth 
is much faster than the diffusion process described above (provided there are particles to be 
absorbed in the neighbourhood of the membrane), and it will stop only when the adsorption 
of particles has effectively changed the chemical potential in the proximity of the membrane. 
We thus obtain the boundary condition of the diffusion problem: 

f^iR) = /^o + (9) 
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An equivalent argument applies for shrinking membranes when jj, < hq. Now it is possible 
to get the chemical potential field for an isolated membrane at an arbitrary distance r from 
its center. Neglecting the thickness of the membrane is equivalent to know the electrostatic 



potential created by a planar metallic disk 



for R < r, where r = {r,6,(f)) in spherical coordinates. From eq. and the rate of 
growth of an isolated membrane of radius R is 



dR _ AaX /TA 1 
H \ Y ~ R 



(11) 



where A = fioo — fio is a measure of the supersaturation of the system. 

We turn now to the study of an ensemble of isolated membranes by introducing a size 
distribution function, fmiR,t) which gives the number of membranes of radius R per unit 
volume at time t. fm{R,t) verifies the continuity equation 

^ + _^ Uy^) = 0, (12) 

dt dR ' ^ ^ 

where vr = dR/dt given by eq. (|lTD . Finally, a closed set of equations is obtained by 
imposing the conservation of the total amount of amphiphiles: 

xA{t) + / dR TTR'TfM t) = Q, (13) 
Jo 

where the first term represents the excess number of particles per unit volume that remain 
in solution {x is the bulk compressibility) and the integral gives the number of particles per 
unit volume that belong to the membranes. With initial conditions for fm{R,t) and A{t) 
the system of equations (0), (|12[) and (|I3|) is fully determined. Following the ideas of the 
LSW theory it is straight-forward to get the asymptotic behaviour of the system (see for 
example |15[). We summarize the main results: 



R{t) = R,{t) = i-^j , (16) 



fMt)-^^P{^y (17) 
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where N{t) is the total number of membranes per unit volume, R{t) is the mean radius and 
P(x) and / are 



P(x) = { (18) 
2 < X, 



oo 

dx x^P{x) = 1.1094 (19) 

It is important to note that within the present model the growth rate R oc t^/^ with a 
growth exponent n = 1/2 instead of 1/3 as predicted by LSW. Such a difference is due to 
the two-dimensional nature of the aggregate, which allows for a faster growth. 
Our asymptotic formulae are expected to be valid for 

/o/ (^*) » 1. (20) 

where Aq is the initial supersaturation of the system. Unfortunately this equation does not 
always hold because before reaching a long enough time some membranes can transform 
into vesicles. Thus, this restriction must hold before Rdt) < 2Rt ~ 8ks/\. From eq. (|16|), 
eliminating t in (|2 



'o.'^ (^) » 1. (21) 

Taking rough estimates of the parameters shown in Table I we get 18 for the left-hand-side 
of eq. (reflog2). Thus, real systems hardly reach this universal distribution function. 



B. Vesicles. 

We now consider the case of a system formed only by spherical vesicles and no planar 
membranes. In the case of an isolated vesicle of radius S, the excess grand potential is: 

An^ = -AnVS'^ifi - /io) + 4:7TKs. (22) 

It is clear from this equation that all vesicles, independently of their size, will grow until the 
chemical potential in the surroundings is /xq- Thus, it is easy to get the three equations that 
define the problem: 

dS _ aA 

~dt~^' ^^^^ 
t + ^(/.^.)-0, (24) 



6 



and 



xA{t) + / dS47rS^rUS,t) = Q, (25) 

<^ 

This set of equations lends itself to an analytical solution in terms of the initial conditions. 
The main result is that the size distribution function does not change its shape, 

f^(S,t) = fo{S-^J^dt'Ait')) (26) 

where /o(>S') represents the initial distribution of vesicles. In this case A{t) decays exponen- 
tially. The reason for this behaviour is clear from eq. (^), which indicates that all vesicles 
grow at the same rate independently of their size. 

It is interesting to notice that eqs. ( PB| - P5D were obtained under the assumption that 
the slowest growth process is the diffusion of particles (diffusion limited growth). But, for 
a vesicle to grow, it is necessary to fill the interior with water, thus a certain amount of 
water needs to overcome an energy barrier when crossing the bilayer structure. Under some 
circumstances this can be the slowest process, and then eq. ( PBj ) should be modified. In this 
case the number of water molecules that cross the bilayer per unit time will be proportional 
to the area of the vesicle. For small enough supersaturations (or large enough curvature 
constant) the shape of the vesicle will remain spherical, and the change in volume is : 

^ = AttS^^ = 47tS^F{A) ^ CS^A (27) 
dt dt ^ ^ ^ ^ 

where F{A) is a function related to the properties of the barrier that verifies F{0) = 
and, thus, it is expected to show a linear behaviour for small supersaturations, C being a 
constant. This equation is equivalent to eq. ( p3D except for a constant, thus the basic result, 
eq. (|26|) still holds except for the time scale. 



C. Membranes and Vesicles. 

When there are no more membranes closing themselves to form vesicles, the distribution 
function does not change in shape, it merely moves to larger sizes until all the excess 
particles in the solution are exhausted. But in order to relate the distribution function when 
the nucleation process finishes it is necessary to allow for the existence of both membranes 
and vesicles and permit the spontaneous transformation from one to the other. We study 
this case assuming that this transformation occurs at one particular radius, Rt = 4:Ks/ A and 
St = Rt/'^, (estimated from the spherical cap model) and the time scale of the process is 
totally negligible compared with diffusion times. Now we have to consider size distribution 
functions for both membranes and vesicles. 

^ + ^ {vRfm) = -S{R - RT)Im (28) 

df d 

+ i^sf.) = S{S - ST)Im (29) 
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POO POO 

xA(t)+ / dRnR'^TfmiR,t) + / dS 47rS^Tf^{S,t) = Q, 
Jo Jo 



(30) 



where Im is the number of membranes per unit volume that transform into vesicles at time 
t, Im = max{0,VRj,fm{RT,t). 

We can change to dimensionless units: 

R = xL, S = xL/2, 
t = tT, 
A = 6D, 

fm = fv = ^F, 

where we have made use of the fact that there are no membranes for R > Rt and there 
are no vesicles for S < St join both size distribution functions into one single function 
ip{x,T): for X < xt = Rt/L it represents membranes of radius R = xL and for x > xt 
vesicles of radius S = xL/2 . With these changes of variable the new equations read: 

i^ + £((^-^H=0 for x<xr 

(32) 

with a special boundary condition at x = xt, and 

6{t) + / dxx^ip{x,T) = 1 (33) 
Jo 

These equations have been solved numerically. We have applied a standard upwind 
algorithm |jT^ discretizing both in time and in space with At = 0.02 and Ax = 0.1. The 
initial conditions were chosen having in mind the particular values of the parameters shown 
in Table I. We consider that initially we only have membranes and assume that the size 
distribution function is gaussian with mean Ri = lOip f» 4 x lO^^m and a width cTj = 
Ri/4 = 10~^m. Our study applies when the nucleation process is already completed. The 
energy barrier for the formation of new membranes must be much lower that /cbT; this means 
that Ao < 7rAV(rfcBT) ^ 8 x 10'^^ J/ molec. We have selected Aq = 3 x 10'^^ J/molec 
which corresponds to a critical radius Rc{0) = 3Ri/A. The only parameter that remains to 
be fixed is the height of the initial distribution, or equivalently the value of Q. We arbitrarily 
consider that at the initial time 90% of the total excess of particles is already in a membrane 
whereas only the remaining 10% is solved in water, leading to Q = 8 x lO^^molec/m^ . The 
results of the calculation are shown in Figure 2. Instead of the size distribution function 
(whose total integral tends to zero in time) we plot x'^ip{x, r) which is proportional to the 
probability of finding a molecule in an aggregate of size x at time r, and its total integral 
tends to one. In the early stage the distribution function moves to larger sizes and spreads, 
decreasing the height of the peak. At r ~ 1300 some membranes start to transform into 
vesicles and the rate increases until r ~ 3000 when it starts to decrease. At r ~ 5200 
the critical radius becomes larger than the transformation radius and ,thus all membranes 
shrink. At that time the final distribution function for vesicles is known; it will just translate 
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(31) 
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to larger sizes absorbing all remaining particles in the solution and in the membranes. In 
Figure 3 is shown the critical radius, Xc versus r. For these initial conditions the behaviour is 
quite monotonous. I is almost linear from r = 500 to r = 6000. For larger times, when most 
of the membranes have disappeared we would reach the expected exponential behaviour. 

A totally different behaviour can be obtained from different initial conditions. The 
results arc shown in Figures 4 and 5. The only difference with respect to the previous 
case is that we have assumed that at the initial time only 5% (instead of 90%) of the total 
excess of particles is forming the initial membranes. In the early time the critical radius is 
lower than the mean radius and there is a large amount of material to be absorbed; this 
permits a very fast growth of the membranes without broadening the distribution function 
and as soon as r ~ 100 there are already some vesicles formed (this early time growth 
process is not shown in the figure because it is covered by later data). This fast growth also 
makes the critical radius increase and, as soon as some membranes are formed, it crosses the 
peak of the distribution of membranes. As the distribution function is highly peaked this 
crossing suddenly makes most of the membranes shrink. For a period of time there is an 
equilibrium where the critical radius remains constant as the vesicles grow at the expense 
of the membranes. When the number of membranes is small the critical radius crosses the 
transformation value. The final distribution function of membranes is highly peaked as a 
result of the very fast early growth. 

III. CONCLUSIONS. 

We have presented a study of the growth processes for amphiphilic membranes. These are 
two dimensional aggregates in a three dimensional bath, with peculiar features: the surface 
vs. volume balance which controls the growth for droplets is changed to a line vs. surface 
balance for planar membranes, or it may be avoided by closed vesicles, without open edge, at 
the price of a size independent curvature energy. Our work here extends the classical Lifshitz- 
Slyozov- Wagner theory to consider this problem. We have shown that, contrary to the 
nucleation of droplets, there is no asymptotic limit in which the size distribution of aggregates 
becomes independent of the initial conditions. As shown by our numerical solutions of the 
coupled equations for planar membranes and spherical vesicles, the asymptotic form of 
the distribution function for vesicles changes dramatically with the initial configuration. 
Allowing for the large uncertainty in the experimental values of several parameters in the 
Table I, and for the overall complexity of the problem, our equations may provide a first 
guide to the systematic understanding of these processes. 

It is important to make a comment regarding the absence of fluctuations in our treat- 
ment. The LSW theory neglects fluctuations because it is expected that in average the 
concentration will follow the gradients in chemical potential created by the existence of 
aggregates. But this assumption conflicts with some other limits taken in the theory in 
particular large times, i.e. large domains, and small supersaturations. In this case the gra- 
dients created by the growth process can be of the same order or lower than the gradients 
created by fluctuations in the concentration. That means that some aggregates could shrink 
even though their radius is larger than the critical radius and this effect can change the size 
distribution function at late enough times. In particular for vesicles, this theory predicts a 
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complete degeneracy in the equilibrium size distribution while fluctuations are expected to 
break this degeneracy imposing the the distribution function which maximizes entropy [18 . 
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TABLES 



TABLE I. Orders of magnitude of some parameters used in the text. Rt has been estimated 
from the spherical cap model, and see text Section II.C for the estimation of Aq. 
keT 4 X 10-21 J 



Pv 


10i« 


molec/rn? 
molec? 1 J 


X 


2.4 X 10^^ 


kg 


10-18 


J 


A 


10-11 


J/m 


r 


1019 


molec/w? 




2 X 10-^ 


m 


Ao 


3 X 10-23 


T/molec 
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FIGURES 



FIG. 1. Curvature energy of a membrane as a function of its shape for a spherical cap model 
(see text). Ti.{A,R) is the curvature energy (in units of Kg) for a membrane of area, A, and edge 

1 /2 

circle of radius, R (in dimensionless units), u = R (tt/A) ' , where u = corresponds to spherical 
vesicles and u = 1 corresponds to planar membranes. The area for the spontaneous shape transition 
is Ac = TT (4«;s/A)^ 

FIG. 2. Time evolution of the size distribution function in dimensionless units, see eq. (|3l|), 
for initial conditions where 90% of the excess particles particles are forming membranes (see text 
for details), x = xt = 66.7 is the radius for the spontaneous shape transition from membrane to 
vesicle. The region x < xt refers to planar membranes and x > xt to spherical vesicles. 

FIG. 3. Critical Radius versus time in dimensionless units, see eq. (|3l|), for initial conditions 
where 90% of the excess particles particles are forming membranes (see text for details). The 
dotted line corresponds to x = xt- 

FIG. 4. Time evolution of the size distribution function in dimensionless units, see eq. (|3l|), 
for initial conditions where 5% of the excess particles particles are forming membranes (see text 
for details), x = xj- = 66.7 is the size of spontaneous shape transition from membrane to vesicle. 
The region x < xt refers to planar membranes and x > xt to spherical vesicles. 

FIG. 5. Critical Radius versus time in dimensionless units, see eq. (|3l|), for initial conditions 
where 5% of the excess particles particles are forming membranes (see text for details). The dotted 
line corresponds to x = xt- 
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